function P = FindQuadrature(im)

chull = regionprops(im, 'ConvexHull');

P = minBoundingBox(chull.ConvexHull');
